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ABSTRACT 

High sensitivity H I observations of the nearby spiral galaxy NGC 2403 ob- 
tained with the VLA are presented and discussed. The properties of the extended, 
differentially rotating H I layer with its H I holes, spiral structure and outer warp 
are described. In addition, these new data reveal the presence of a faint, extended 
and kinematically anomalous component. This shows up in the H I line profiles as 
extended wings of emission towards the systemic velocity. In the central regions 
these wings are very broad (up to 150 km s^^) and indicate large deviations from 
circular motion. We have separated the anomalous gas component from the cold 
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disk and have obtained for it a separate velocity field and a separate rotation 
curve. The mass of the anomalous component is 1/10 of the total H I mass. 
The rotation velocity of the anomalous gas is 25—50 km s~^ lower than that of 
the disk. Its velocity field has non-orthogonal major and minor axes that we 
interpret as due to an overall inflow motion of 10—20 km s~^ towards the centre 
of the galaxy. 

The picture emerging from these observations is that of a cold H I disk sur- 
rounded by a thick and clumpy H I layer characterized by slower rotation and 
inflow motion towards the center. The origin of this anomalous gas layer is un- 
clear. It is likely, however, that it is related to the high rate of star formation in 
the disk of NGC 2403 and that its kinematics is the result of a galactic fountain 
type of mechanism. We suggest that these anomalous H I complexes may be 
analogous to a part of the High Velocity Clouds of our Galaxy. 

Subject headings: galaxies: individual (NGC 2403) — galaxies: structure — 
galaxies: kinematics and dynamics — galaxies: halos — galaxies: ISM 

1. Introduction 

In recent years the study of the vertical structure and kinematics of the H I disks of 
spiral galaxies has been pursued almost exclusively with observations of galaxies viewed 
either edge-on or face-on. The study of edge-on systems has revealed vertical extensions of 
the H I layer up to several kpc from the plane and, in the case of NGC 891, also a slower 
rotation velocity of the gas in the halo (Swaters, Sancisi, & Van der Hulst 1997). In the 
face-on or nearly face-on galaxies the observations have shown vertical motions of neutral 
gas frequently associated with 'holes' in the H I distribution (Puche et al. 1992; Kamphuis 
1993). These results point at a substantial gas circulation between disk and halo (review by 
Sancisi (1999)). The ejection of gas out of the disk is thought to be related to star formation 
in the disk and to be due to stellar winds and supernova explosions that blow out ionized 
gas. After coohng, this gas falls down to the plane. Such processes are generally described 
as galactic fountains (Shapiro & Field 1976; Bregman 1980). There are cases, however, in 
which the vertical motions originate from external events such as collisions with intergalactic 
clouds and with small companions (van der Hulst & Sancisi 1988). Accretion of intergalactic 
'primordial' gas, as proposed for the High Velocity Clouds in our Galaxy (Oort 1970), may 
also play a role. 

The present study represents a flrst attempt to use objects at intermediate inclination 
angles which offer the advantage that, for a given object, information is obtained on both 
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the density structure and the kinematics of H I in the vertical direction. The nearby spiral 
galaxy NGC 2403 is a suitable candidate because of its inclination angle of ~ 60°, and also 
because of its extended H I layer (about twice the optical size), its regular kinematics and 
its symmetric, flat rotation curve. NGC 2403 is an isolated ScIII galaxy (Figure 1) with 
a bright optical disk of about 8 kpc diameter surrounded by outer, faint-surface-brightness 
'fossil arms' (Sandage, Tammann, & Van den Bergh 1981). It hosts rich H II regions with 
brightnesses comparable to that of the 30 Doradus complex (Drissen et al. 1999). 

Previous H I observations with the Westerbork Synthesis Radio Telescope (WSRT) 
(Begeman 1987; Sicking 1997; Schaap, Sancisi, & Swaters 2000) of NGC 2403 showed sys- 
tematic asymmetries in the H I velocity proflles along the major axis in the form of wings 
on the side of the lower rotation velocities, towards systemic. This pattern (we refer to it as 
'the beard') is similar to that found in galaxies viewed edge-on or observed with low angular 
resolution but it is totally unexpected in NGC 2403 if its H I disk is thin. Schaap et al. 
(2000) re-analyzed Sicking's (1997) observations of NGC 2403 and investigated the effect of 
the layer thickness in order to explain such a pattern. They concluded that the beard may 
be produced by H I gas located above the plane (in a thick disk or halo) and rotating more 
slowly than the disk. Their picture is very similar to that proposed by Swaters et al. (1997) 
for NGC 891. 

We obtained new H I observations of NGC 2403 with the Very Large Array.^ These 
have led to a considerable improvement of the observational picture and to the detection of 
new, very faint emission from H I with remarkably anomalous kinematics. 

2. Observations 

We have observed NGC 2403 with the CS configuration of the VLA on four separate days 
in January and February, 1999. This configuration differs from the regular C-array in that 
the middle (fifth) telescope on the North arm is located at the innermost North arm location 
of the D array. This leads to a spatial resolution similar to that of the regular C-array, while 
adding short baselines - and hence sensitivity at larger angular scales - not-found in regular 
C-array observations (Rupen 1998). Each observation was 12 hours in duration, of which 
approximately 10 hours were spent on source. The remaining time was used to observe one 
phase cahbrator (0841+708) and three flux calibrators (3C48, 3C147, and 3C28). The latter 
were used for bandpass calibration as well. The observational parameters are summarized 
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in Table 1. 



2.1. Data reduction 

Calibration and editing of the data were done separately for each of the four data-sets, 
using standard procedures of the AIPS (Astronomical Image Processing System) package. 
Nineteen channels at both ends of the band (excluding the very edges) were identified as 
being free of line emission. The AIPS task UVLIN (Cornwell, Uson, & Haddad 1992) was 
then used to interpolate this continuum emission across the range of channels containing line 
emission, and subtract it from the uv data. 

The resulting continuum-free uv dataset containing the combined data for all four days 
was then further reduced using the MIRIAD (Multichannel Image Reconstruction, Image 

Analysis and Display) package. The data cube was made using Robust or Briggs' weighting 
(Briggs 1995) with robustness 0.2 getting a spatial resolution of ~ 15" or 230 pc at the 
adopted distance of 3.18 Mpc (Madore & Friedman 1991). Another data cube was made by 
applying a Gaussian taper in the visibility domain. This data cube has a spatial resolution 
of ~ 30". The two cubes have 62 channels ranging in velocity from —24.7 km s^^ to 289.6 
km s~^ and spaced by 5.15 km s~^. To these data a Manning smoothing was applied in 
velocity, leading to a velocity resolution of 10.3 km s~^. The dirty images were deconvolved 
using a Clark CLEAN algorithm (Clark 1980). To identify the region of emission, the cubes 
were firstly CLEANed with a cutoff of 3 times the r.m.s. noise. The resulting images were 
then used to interactively define the regions of line emission in each channel. Next, a deep 
CLEAN (cutoff of ~ 0.2 a) was performed. 

After cleaning, a few channel maps still showed some spurious effects (Figure 2, top left 
corner of maps from 217.5 km s~^ to 248.4 km s~^). These may indicate the presence of very 

extended H I associated with the outer parts of the warp. Some channel maps show serious 
contamination from H I emission from our Galaxy at velocities around km s^^ (Figure 2, 
see the channel maps between ~ —10 km s~^ and ~ +10 km s~^). This contamination is 
slightly reduced by applying the Banning smoothing and also by using high resolution data 
because it generally has a greater angular scale than the emission from NGC 2403. 

Comparison of our global B I profile with a previous single dish profile from Rots (1980) 
(Figure 1) shows we are not missing significant B I fiux due to missing short spacings. The 
final data cube at ~ 15" resolution has a r.m.s. noise per channel of 0.17 mJy beam~^, 
while the lower resolution data cube has a r.m.s. noise of 0.22 mJy beam~^. The minimum 
detectable column density (5 a detection in one velocity resolution element in the 30" reso- 
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lution cube) is 2.0 x 10^^ cm~^ (0.15 M© pc~^) corresponding to ~ 2.6 x 10^ M© beam~^. 
This very high sensitivity has permitted us to detect very faint components in the neutral 
hydrogen that were not seen in previous observations. A summary of the parameters for the 
two cubes is given in Table 2. 



2.2. Data description 

Once the data cubes were obtained, the subsequent analysis was done using the GIPSY 
(Groningen Image Processing SYstem) package. Figure 2 shows every second channel map 
of the 30" resolution cube of NGC 2403. We present here the smoothed data because they 
show more clearly the faint H I components. Some representative channels at high (15") 
resolution arc shown in Figure 8 overlaid on an optical DSS image. Wc also present several 
position-velocity (p-v) diagrams parallel to the major (Figure 4) and minor (Figure 5) axes. 
The cuts are shown in Figure 3 overlaid on the total H 1 map. 

The full resolution total H I map (Figure 1) was obtained by adding the channel maps 
containing H I emission (from —19.6 km s~^ to 279.3 km s~^). For each channel map the 
areas with H I emission were outlined using masks made on the 30" resolution smoothed 
data cube. The total H I mass, corrected for the primary beam attenuation, is 3.24 ± 0.05 
X 10^ Mq. The value of the mass derived for this galaxy by Begeman (1987) was 3.0 x 10^ 
Mq while Sicking (1997) with more sensitive observations derived 3.13 x 10^ Mq assuming 
the distance of 3.18 Mpc. For comparison the single dish mass (Rots 1980) is 3.31 x 10^ 
Me- 

To first order, the distribution and the kinematics of the H I are regular. Figure 6 (left 
panel) shows the total H I at low (60") resolution overlaid on the DSS optical image. The 
radius of the H I disk down to a column density of ~ 0.2 MqPc^^ (2.5 x 10^^ cm~^) is 
about 22 kpc (the Holmberg radius is 13 kpc). The H I disk appears to be full of holes, some 
having kpc sizes (Mashchenko, Thilker, & Braun 1999). The very outer layer is warped. The 
warp is highly asymmetric, more pronounced on the eastern side of the galaxy (see channels 
from 176.2 km s^^ to 227.8 km s~^) than on the western side where it is just faintly visible 
(channels from 42.3 km to 62.9 km s~^). The warp is also visible in the total H I map and 
in the velocity field at low resolution shown in Figure 6 where the H I disk on the N-E side 
appears to be more extended. A summary of the optical and H I parameters for NGC 2403 
is given in Table 3. 
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2.2.1. The anomalous gas 

The most striking characteristics of the kinematics of the neutral hydrogen in NGC 2403 
are illustrated in Figure 7. This diagram offers a display of the H I emission along the major 
axis of NGC 2403. It was obtained by integrating the 15" full resolution data over a slice 1' 
wide in order to improve the signal to noise ratio. Clearly, the H I line profiles are strongly 
asymmetric with respect to the peak of the line emission. Note that the white squares, which 
mark the rotation curve, follow closely the ridge of the H I. There are H I tails extending 
(as a beard) towards the systemic velocity. The brighter of these tails had been seen in 
previous 21-cm observations with the WSRT and already reported by Schaap et al. (2000). 
Because of the higher sensitivity of the present observations the tails are now much easier to 
recognize and they appear to extend over a much larger region of the position-velocity map. 
A clearer picture of the extent, the structure and the kinematics of the gas in the beard is 
now emerging. The density distribution is characterized by sub-structures in the form of 
filaments and spurs. These are observed all over the H I disk at velocities that differ tens of 
km s~^ from rotation. The most remarkable filament has a size of ~ 8 kpc and is visible in 
the channel maps between 93.8 km and 155.6 km s~^. It has a coherent structure and 
kinematics and its total H 1 mass is ~ 1 x 10^ Mq. It does not show any obvious relation 
with any feature in the optical disk (see upper right and bottom left panels of Figure 8) 
or with the H I spiral arms. The —V and —2' p-v plots parallel to the major axis (and 
partly also that along the major axis itself) in Figure 4 show this filament elongated on 
the low velocity side between 0' and —8'. Remarkably, it follows the overall pattern of the 
low-density gas suggesting that it is just a prominent condensation of gas within it and not 
a separate feature. 

The p-v plots parallel to the minor axis (Figure 5) show an asymmetric distribution of 
the low-level emission (tracing the beard). In the diagrams from —4' to —1', at velocities 
below the systemic velocity (horizontal line at 133 km s"^), this low-level emission appears 
preferentially in the left part of the diagrams, between 0' and ~ -|-10' (North-East side of 
the galaxy), whereas at velocities above systemic it is more visible in the right part of the 
plots. In the diagrams from +1' to +4', at velocities above systemic the low-level emission 
is more prominent on the right part of the diagrams; at velocities below systemic it is in 
the left part. A similar asymmetric distribution of the low-level H I emission is also seen in 
the channel maps close to the systemic velocity (Figure 2) where the lower contours do not 
follow the H I ridge and seem displaced counter-clockwise. This effect indicates a rotation 
of the kinematical axes of the faint emission with respect to the axes of the cold disk. As we 
will discuss later (see Section 3.2.1) such a change in the kinematical position angles can be 
attributed to an overall radial motion of the anomalous gas of about 10—20 km s~^ towards 
the centre of the galaxy. 
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Apart from the extended low-level emission described above, the p-v diagrams in Fig- 
ure 4 and Figure 7 also show emission in the forbidden quadrants. This is particularly clear 
in the upper right quadrant in Figure 7 (0' to —4' North- West with respect to the centre) 
along the major axis at projected velocities (around 200 km s~^) which differ from rotation 
by more than 130 km s~^. The spatial location of this gas can be seen in the channel maps 
between 165.9 km s"^ and 217.5 km s~^ where it is visible just North- West of the centre. 
The total mass of this gas is 4.3 x 10^ Mq. The channel maps between 62.9 km s~^ and 
104.1 km s"^ show a faint counterpart on the South-East side of the galaxy with a total 
mass of 1.7 x 10^ Mq. All this gas at extreme anomalous velocities (the 'forbidden' gas) 
seems to be confined to the region corresponding to the bright optical part of NGC 2403 (see 
Figure 8, top left and bottom right panels). Its overall kinematical pattern, as it appears 
in Figure 7, indicates that all the anomalous gas, including the forbidden one, forms one 
coherent structure, which closely follows the rotation of the system. 

In the following we refer to all the gas at projected anomalous velocities with respect 
to the rotation curve (beard -|- forbidden gas) as the 'anomalous gas'. 

3. The separation of the anomalous gas 

We have attempted to separate the anomalous from the regular H I assuming that the 
H I profiles are formed by two components: 1) a narrow one of Gaussian shape centred on 
the rotation velocity of the galaxy (the cold disk), and 2) a broader one, mainly at lower 
rotation velocities and with unknown profile shape (the anomalous gas). We have modelled 
the cold disk by fitting a Gaussian function to the line profile after clipping the data at 30% 
of the peak. With this clipping we fit only the upper part of the profiles. The differences 
between the Gaussian fit with and without the clipping are about 1 km s~^ for the value 
of the rotation velocity and about 3 km s~^ for the velocity dispersion. The velocity field 
and the distribution of the velocity dispersion of the cold disk (Figure 1), obtained with this 
Gaussian fit, have then been used to generate a 'Gaussian cube' that has been subtracted 
from the data leaving the residual anomalous gas. A second method has been tried to 
generate the Gaussian cube for the cold disk by fitting half a Gaussian to the high velocity 
side of the line profiles. This is based on the assumption that the anomalous gas does not 
contribute to the high velocity side of the line profile. The obtained half-Gaussian has then 
been folded on the low velocity side of the profile and subtracted from the data. The results 
obtained with the two methods are very similar. 
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3.1. The cold disk 

The velocity field and the map of the velocity dispersion of the cold disk are shown in 
Figure 1. The velocity field is regular and symmetric. The velocity dispersion has average 
values around 9—10 km s~^ except in the central region and in the spiral arms where it reaches 
values of up to 10—15 km s"^ (darker grayscale). After correction for the instrumental 
resolution these values are about 1 km s~^ lower. In order to determine the kinematical 
parameters and the rotation curve of NGC 2403 we have performed the standard tilted ring 
fit of the velocity field (Begeman 1987). The rings have been fitted with a radial increment 
of 15". Points were weighted by the cosine of the azimuthal angle with respect to the major 
axis. In the determination of the central position and the systemic velocity only the rings 
with radii smaller than 13 kpc were considered. The values derived for the position of the 
centre and for the systemic velocity are reported in Table 3. While keeping these parameters 
fixed, we proceeded to fit the position and inclination angles. The position angle (Figure 9) 
shows a clear deviation from the mean value in the central region between 2 and 4 kpc and 
remains constant out to about 15 kpc where it starts decreasing in the region of the warp. 
The inclination angle shows variations in the central region around a mean value of about 
61°. Then it rises to more than 65° beyond 15 kpc. The mean (error weighted) values for 
the position and the inclination angles were derived between 4 and 15 kpc and resulted in 
124.5 ± 0.6° and 62.9 ± 2.1° respectively. The values obtained by Sicking (1997) were 123.9 
and 60.9 degrees respectively while Begeman (1987) found 122.5 and 60.2. The lower values 
for the inclination angle found by these authors agree with the value we get considering only 
the inner regions (radii < 10 kpc) of the disk. 

While keeping the other parameters fixed, the rotation curve was derived separately for 
the approaching and receding side of the galaxy, as shown in the top left panel of Figure 9. 
The line shows the rotation curve obtained by Begeman (1987). Points beyond 20 kpc were 
determined using the velocity field at 30" resolution. Because of the large deviations of 
position and inclination angles from their mean values beyond 20 kpc they were left free 
to vary in the fit. The rotation curves for the two sides indicate a good overall symmetry 
except for the regions beyond 16 kpc where the approaching side has velocities 5—10 km s"^ 
higher than the receding one. 

Figure 9 also shows the radial column density fitted with an exponential law for radii 
larger than 4 kpc (the inferred H I scalelength is 5.7 ± 0.1 kpc), and the velocity dispersion 

derived with the formula adisp^ = o'ohs'^ + o"instr^ where Uinstr is the instrumental velocity 
resolution (~ 2 x channel separation) and (Jobs is the a of the Gaussian function fitted to 
the clipped line profile. The derived velocity dispersion shows a decrease from the central 
values of ~ 12 km s^^ to values of ~ 8 km s~^ in the outer regions. 
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3.1.1. The mass model 

The rotation curve of NGC 2403 was derived again by applying the tihed ring fitting to 
the velocity field of the whole galaxy. For points beyond 20 kpc we again used the smoothed 
data cube with no fixed position and inclination angles. As an estimate of the errors we took 
the difference in velocity between the approaching and receding side. We have subsequently 
performed the standard analysis of the rotation curve with the decomposition into three 
mass components (Begeman 1989): the stellar disk, the gas disk, and the dark matter halo. 
No stellar bulge is visible in the optical luminosity profile (Kent 1987). For the stellar disk 
we have taken the brightness profile published by Kent (1987) while for the H I profile we 
multiplied our data by a scale factor of 1.4 to take into account the Helium abundance. We 
assumed a value of 0.4 kpc for the scale height of the stellar layer (van der Kruit & Searle 
1981) and 0.2 kpc for gas layer (Sicking 1997). The dark matter halo was modelled with 
the spherically symmetric density distribution p{R) = po{l + -^)~'^, where po is the central 
density and Rq is the core radius. 

In Figure 10 we show four different fits. The maximum disk model gives a scaling factor 
(M/L) for the stellar component of 2.3, the 'Bottema disk' model (Bottema 1993) has a 
fixed M/L of 1.4, and the third model is the so-called maximum halo (no stellar disk). The 
fourth panel shows a model with no dark matter halo and a free scale factor for stars and 
H I; we obtained an M/L of 2.3 and a value of 9.9 for the H I scahng factor, similar to that 
found by Hoekstra, Van Albada, & Sancisi (2001). The first three fits are comparable. Note 
that the last model gives a higher discrepancy in the outer regions than found by Hoekstra et 
al. 2001. However, this latter produces the best fit of the inner substructures of the rotation 
curve. The present rotation curve is more extended (about 2 kpc) than that by Begeman 
(1987) and does not show any sign of declining. The total dynamical mass inferred within 
a galactic radius of 22.5 kpc is 9.4 ± 0.7 x 10^° Mq. This value is consistent with that 
obtained by Begeman (1987) (8.4 ± 0.4 x 10^° M©) for a galactic radius of 20 kpc. 



3.2. The anomalous gas 

Figure 11 shows the p-v diagram for the anomalous gas along the major axis of NGC 2403 
after the subtraction of the cold disk. The filled squares show the rotation curve derived for 
the cold disk. Figure 12 shows the distribution of the anomalous gas overlaid on the optical 
picture (left panel). Its outer radius is about 15—16 kpc and the total H I mass is ~ 3 x 
10^ Mq (1/10 of the total H I mass of NGC 2403, 0.3% of the total dynamical mass). For 
the velocity dispersion of the anomalous gas the data indicate values of 20—50 km s~^. 
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Where is this anomalous gas located in space with respect to the cold disk? Is it 
located above the disk, in the halo region? The vertical extent of this gas layer is difficult to 
determine precisely because the inclination of the galaxy is not sufficiently high. Schaap et 
al. (2000) using a 3D modelling of the gas layer arrived at possible values for its thickness 
between 1 and 3 kpc (FWHM). Similar values were found for the edge-on galaxy NGC 891 
(Swaters et al. 1997) whose H I 'halo' is probably a component similar to the anomalous gas 
in NGC 2403. A similar modeUing performed by us gives an upper limit for the thickness of 
the anomalous gas of about 3 kpc, while the thickness of the disk is presumably lower than 
0.4 kpc (Sicking 1997). Also from a dynamical point of view this separation into two systems 
of different thicknesses is natural as the anomalous gas and the cold disk, being coUisional 
systems with separate kinematics, can not spatially co-exist in the same thin layer. 



3.2.1. Velocity field and radial motions 

Figure 12 (right panel) shows the intensity- weighted velocity field of the anomalous 
gas. Clearly the kinematics is dominated, as for the cold disk, by differential rotation. The 
projected kinematical axes, however, appear to be rotated counter-clockwise with respect to 
those of the cold disk. This is more evident for the minor axis (thick line). Furthermore, 
the two axes are not orthogonal. This is a result of the asymmetric distribution of the 
anomalous gas already noted in the description of the p-v parallel to the minor axis of the 
galaxy (Section 2.2.1). A simple explanation is that of an overall inflow of the anomalous 
gas towards the center of the galaxy (Fraternali et al. 2001). 

We have explored the possibility of such a radial flow quantitatively using a tilted ring 
fit to the velocity field of the anomalous gas. This was done by keeping the centre of the 
galaxy and the systemic velocity fixed to the values found for the cold disk. The resulting 
inclination angle is similar to that of the disk, but the position angle, as expected from the 
counter-clockwise rotation of kinematical axes, turns out to be significantly higher (~ 131°). 
Assuming such a rotation of the kinematical axes is indeed due to radial motions, we have 
included in the tilted ring model a radial component for the gas velocity. In order to decide 
about the sign of such radial motion —inflow or outflow— it is necessary to know which is the 
near and which is the far side of the galaxy. We have taken the South- West side as the near 
side of NGC 2403. This is based on the assumption that the spiral arms are trailing and is 
also supported by the presence of dust absorption features (sec optical images in Sandagc et 
al. 1981). Keeping the position angle fixed to the value found for the cold disk (124.5), the 
fit gives a radial flow velocity of —13.6 ± 1.2 km s~^ (the minus sign means inflow) from 4 
to 15 kpc. 
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ls it possible that also the cold disk has a non-zero radial velocity? We tested this 
possibility by considering both radial velocity and position angle as free parameters in the 
tilted ring fit to the velocity field of the cold disk. Figure 13 (right panel) shows the trend 
of radial velocity for both the disk and the anomalous gas obtained with this method. The 
mean values of radial velocity between 4 and 15 kpc are —3.1 ± 0.3 km s^^ and —16.2 ± 1.1 
km s~^ for the disk and the anomalous gas respectively. Thus the velocity field of the cold 
disk also shows a small distortion that could be explained as a weak radial inflow. However, 
an alternative and more likely interpretation for this distortion is that of triaxiality of the 
potential. Indeed, the deviation from circular motion found here is comparable to that found 
by Schoenmakers, Pranx, & de Zeeuw (1997) (about 2 km s~^) with a Fourier analysis of 
the orbits. Such a weak triaxiality would, instead, be totally insufficient to account for the 
effect (a factor of 10 greater) observed for the anomalous gas. 

Figure 13 also shows the rotation curve of the disk compared with the mean rotation 
velocity of the anomalous gas. The difference is about 25 km in the outer parts, increasing 
to 40—50 km s~^ in the inner regions. Table 4 summarizes the properties of the H I disk and 
of the anomalous gas for NGC 2403. 

4. Summary and concluding remarks 

We have presented new, deep H I observations of the spiral galaxy NGC 2403. These 
have been used to study the structure of the gas layer and the dynamics of the system. The 
main results are briefly summarized here. 

A new faint, extended gas component with anomalous velocities has been detected. This 
H I component (the 'anomalous gas') shows up in the fine proflles as broad wings at velocities 
lower than rotation and partly also at forbidden (non-circular) velocities. Such wings are not 
explained by inclination effects, pure thickness or poor angular resolution, but result from 
the presence of a thick layer of H I which rotates more slowly than the gas in the plane (see 
also Schaap et al. 2000). 

We have separated the anomalous gas from the thin disk of cold H I assuming a Gaussian 
profile for the latter. The anomalous H I has a mass of about 3 X 10^ Mq (1/10 of the total 
H I mass), a mean rotation velocity 25—50 km s~^ lower than that of the disk, and a radial 
infiow of 10-20 km s-\ 

The origin of this gas is still a matter of debate. The main issue is whether it is the result 
of processes taking place in the system itself or whether it could be infall of extragalactic, 
possibly primordial gas (cf. Oort 1970). The overall kinematical pattern of the anomalous 
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gas, its highly coherent structure and its close connection with the disk of NGC 2403 seem to 
point at an internal origin. Indeed, there may be a close relation with the process of active 
star formation going on in the disk of NGC 2403 and a galactic fountain type of mechanism 
may provide the explanation for both its origin and its dynamics. However, the detection 
of H I at forbidden velocities (more than 130 km s~^ from rotation) and the presence of 
coherent structures with sizes up to 8 kpc, are puzzling. Considering the large velocity 
dispersion (about 20—50 km s~^) and the lower rotation velocities found for the anomalous 
H I, it is tempting to describe what we see here in terms of a large asymmetric drift of a 
thick and clumpy gas layer. 

The results of the present H I observations of NGC 2403 have also interesting impli- 
cations regarding the High Velocity Clouds observed in the Milky Way (Wakker & Van 
Woerden 1997) whose distances, nature and origin are still a matter of debate. It is possible 
that the anomalous H I complexes discovered in NGC 2403 are analogous to the HVCs seen 
in the Galaxy, or at least to a class of them. Consider for instance the 8 kpc filament, the 
largest structure found in NGC 2403, with a projected velocity difference from the rotation 
of 60—100 km s^^ and an H 1 mass of about 1 X 10^ Mq. If viewed from inside NGC 2403 
it would probably appear to have a large angular size similar to that of some extended HVC 
complexes e.g. the well known complex C (Wakker et al. 1999). To be as massive as the 
filament in NGC 2403, complex C would have to be as far as 15 kpc from us, compatible 
with the observational lower hmits (Wakker & Van Woerden 1997). Such a distance would 
also be of the same order as the distance (between 4 and 10 kpc) estimated for complex 
A (van Woerden et al. 1999) that shows several similarities with complex C. The fact that 
NGC 2403 is a normal, non-interacting spiral galaxy suggests that the H I features discovered 
here may be common among spiral galaxies and that, perhaps, they have not been detected 
yet because of the low sensitivity of previous surveys. 

Apart from the investigation of the newly discovered anomalous gas component, the 
present observations have also been used for the study of the dynamics of NGC 2403. The 
main results are the finding of an extended, highly asymmetrical warp of the H I outer 
layer and the derivation of a new rotation curve. This has been analyzed with the standard 
decomposition technique and the main results on the disk and the dark matter halo are 
similar to those already obtained by Begeman (1987). 

We acknowledge financial support from the Italian Ministry for the University and 
Scientific Research (MURST). 
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Fig. 1. — From upper left: optical image (DSS) of NGC 2403, global H I profile (dots) 
compared with previous single dish observations (line) by Rots (1980), total H I map, velocity 
field, radio continuum, and velocity dispersion field. These are all at full resolution (15"). 
The column density in the total H I map ranges from 1x10^° to 5x10^^ cm~^. In the velocity 
field (central right panel) the contours are separated by 20 km and the thicker line marks 
the systemic velocity (133 km s~^). The receding side is darker. 
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Right Ascension (J2000) 

Fig. 2. — H I channel maps at 30" resolution (Manning smoothed). The hehocentric radial 
velocities (km s~^) are shown in the upper left corners. The contours are —0.5, 0.5, 1, 2, 5, 
10, 20, 50 mJy/beam; the r.m.s. noise is 0.22 mJy/beam. The cross indicates the kinematical 
centre of the galaxy. 
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Figure 2.2 (continue) 
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Figure 2.2 (continue) 
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Figure 2.2 (continue) 



-20- 




- 21 - 



65''28 




Right Ascension (J2000) 



Fig. 3. — Total H I map showing the cuts parallel to the major and the minor axes used for 
the position- velocity maps displayed in Figures 4 and 5. 
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Fig. 4. — Position-velocity diagrams (30" resolution, Hanning smoothed) parallel to the 
major axis of NGC 2403 (p. a. =124. 5°, see Figure 3). The central panel is the p-v plot along 
the major axis itself, the +1', +2' are from cuts taken North-East respectively at 1' and 2' 
from the centre of the galaxy, the —1', —2' are located South- West (1' is ~ 1 kpc). The left 
side of the diagram corresponds to the South-East side of the galaxy. The central horizontal 
line shows the systemic velocity, the two lower lines mark the channel maps contaminated 
by emission from the Milky Way. Contours are: —0.45, 0.45, 1, 2, 4.5, 10, 20, 45 mJy/beam. 
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Fig. 5. — Position-velocity diagrams (30" resolution, Hanning smoothed) parallel to the 
minor axis of NGC 2403 (p. a. =34. 5°, see Figure 3). The central panel is the p-v plot along 
the minor axis, the +1', +2'... are from cuts taken North-West respectively at 1' and 2'... 
from the centre of the galaxy, the —1', —2'... are located South-East. The left side of the 
diagram corresponds to the North-East side of the galaxy. The central horizontal line shows 
the systemic velocity, the other two horizontal lines mark the channels contaminated by H I 
emission from the Milky Way. Contours are: —0.45, 0.45, 1, 2, 4.5, 10, 20, 45 mJy/beam. 
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Fig. 6. — Total H I overlaid on the DSS optical image and velocity field for NGC 2403 at 
low (60") angular resolution. The presence of the outer warp is visible on the N-E side of 
the disk. The contours are 2.4, 6.4, 13, 26, 52, 78, 126, 207 xlO^^ cm^^ (the lower level is 
2 in terms of r.m.s. noise). In the velocity field the contours are separated by 30 km s~^ 
and the thicker one marks the systemic velocity (133 km s~^). The receding side is darker. 
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Fig. 7. — Position velocity diagram (slice ~ 1' wide) along the major axis (p. a. =124. 5°) of 
NGC 2403. The spatial resolution is 15", the velocity resolution 10.3 km s~^. The central 
horizontal line shows the systemic velocity, the other two horizontal lines mark the channels 
contaminated by H 1 emission from the Milky Way. Contours are: —0.26, 0.26, 0.5, 1, 2, 
5, 10, 20 mJy/beam. The r.m.s. noise is 0.17 mJy/beam. White squares mark the rotation 
curves for the two sides of the galaxy. 
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Fig. 8. — Four representative channel maps at 15" resolution overlapping the optical image 
of NGC 2403. Contour levels are: —0.5, 0.5, 1, 2, 5, 10, 20 mJy/beam. The r.m.s. noise 
is 0.17 mJy/beam. The lower right panel shows the forbidden gas in the North- West side 
of the galaxy while the upper left panel shows the faint emission of forbidden gas in the 
South-East side. The other panels show the bright 8 kpc filament in the West side of the 
H I disk. Note its location outside the bright part of the optical disk. 
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Fig. 9. — Radial profiles of parameters of the tilted ring model for NGC 2403. In the top 
left panel the rotation curves obtained for the two sides of the galaxy are shown. The line 
shows the rotation curve obtained by Begeman (1987). 
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Fig. 10. — Four different mass models for NGC 2403. The filled squares show the observed 

rotation curve, the contributions by gas (dashed), stars (short dashed), and DM halo (long- 
short dashed) and the total rotation velocity (solid). Values of and scaling factor (M/L) 
of the stellar profile are given as well. 
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Fig. 11. — Position-velocity diagram for tlie anomalous gas along the major axis of 
NGC 2403. Contours are 2, 4, 8, 20, 40 in units of r.m.s. noise. The filled squares mark the 
rotation curve of the cold disk. 
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Fig. 12. — Total H I map and weighted mean velocity field of the anomalous gas in NGC 2403. 
In the left panel the density distribution of the anomalous gas is overlaid on the optical DSS 
image of the galaxy. The contours are 2, 6, 15 in units of r.m.s. noise. In the velocity 
field (right) the contours are separated by 30 km s~^ and the thick one shows the systemic 
velocity (kinematical minor axis). Note that this latter is rotated counter-clockwise (by 
about 20 degrees) with respect to the minor axis of the velocity field of the cold disk shown 
in Figure 1. 
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Fig. 13. — Rotation curves (left) and radial curves (right) for the the cold disk (filled squares) 
and for the anomalous gas (open squares) . The mean values of the radial infiow velocity for 
the anomalous gas and for the disk (fines) are respectively 16.2 km s~^ and 3.1 km s~^. 
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Table 1. Observational parameters for NGC 2403. 



NGC 2403 



Observation dates 


Jan 1, 28, 31; Feb : 


Length of observation (hours) 


48 


Time on source (hours) 


40 


Configuration 


CS 


Number of antennas 


27 


Pointing R.A. (J2000) 


07 36 54.5 


Pointing Dec (J2000) 


+65 35 20.0 


Central velocity (km s^^) 


130 


Central frequency (MHz) 


1419.747 


Total bandwidth (MHz) 


3.125 


Total bandwidth (km s~^) 


660 


Number of channels 


127 


Channel separation (kHz) 


24.4 


Channel separation (km s~^) 


5.1 
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Table 2. Parameters of the data cubes. 





Pull resolution 


Smoothed 


Used value of ROBUST 


0.2 





Used taper (") 


no taper 


27 


HPBW (") 


15.2 X 13.3 


29.7 X 29.3 


P. A. of synthesized beam (°) 


91.7 


98.7 


Beam size (pc) 


234 X 205 


458 X 452 


Number of channels 


62 


62 


R.m.s. noise per channel: 






No Manning smoothing (mjy/beam) 


0.28 


0.35 


Hanning smoothing (mJy/beam) 


0.17 


0.22 


Minimum detectable column density: 






(cm-2) 


5.0 X 10^^ 


2.0 X 10" 


(M0/pc2) 


0.38 


0.15 


Minimum detectable mass (Mq /beam) 


2.1 X 10* 


2.6 xlO-* 


Conversion factor mJy/K 


3.0 


0.7 
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Table 3. Optical and radio parameters for NGC 2403. 





NGC 2403 


ref 


Morphological type 


Sc(s) III 


1 


Optical centre (a, <5 J2000) 


7 36 51.92 65 36 0.6 


2 


Kinematical centre {a, d J2000) 


7 36 50.66 ± 1.48 65 36 2.2 ± 5.5 


3 


Distance (Mpc) 


3.18 


4 




7.9 ± 0.7 X 10^ 


5 


Disk scale length (kpc) 


2.0 


6 


R25 (kpc) 


8.2 


5 


Holniberg radius (kpc) 


13.4 


7 


H I scale length (kpc) 


5.7 ± 0.1 


3 


H I disk radius (kpc) 


22.5 


3 


Systemic velocity (km s~^) 


133.2 ± 2.2 


3 


H I total mass (Mq) 


3.24 ± 0.05 X 10^ 


3 


Mean H I inclination (°) 


62.9 ± 2.1 


3 


Mean P. A. (no warp) (°) 


124.5 ± 0.6 


3 


Total mass (Mq) 


9.5 ± 0.7 X 10^° 


3 



References. — (1) Sandage, Tammann, & Van den Bergh (1981) (2) Gallouet, 
Heidmann, & Dampierre (1973); (3) this work; (4) Madore & Friedman (1991); 
(5) de Vaucouleurs, de Vaucouleurs, & Corwin (1976); (6) Wevers, Van der Kruit, 
& Allen (1986); (7) Holmberg (1958). 
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Table 4. Comparison between disk and anomalous gas. 





Cold disk 


Anomalous gas 


Mass (Mq) 


3 X 10^ 


3 X 10* 


Diameter (kpc) 


46 


32 


Thickness (kpc) 


0.4 


<3 


Maximun rotation velocity (km s^^) 


130 


110 


Mean ratlial inflow velocity (km s^^) 


0-3 


10-20 


Mtvm \'('iocil>" (iisixTsiiin (km s ) 


S-12 


20-.^)l) 



